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Abstract 

A  geometrical  theory  is  developed  for  the  analysis  of  surface  wave 
excitation  and  propagation.  The  surfaces  along  which  the  surface  waves 
propagate  may  be  either  curved  or  flat,  and  may  have  either  constant  or 
variable  properties.  The  theory  is  based  on  the  concept  of  a  complex  or 
imaginary  ray.   The  excitation  coefficient  which  enters  the  theory  is  deter- 
mined from  the  solution  of  a  canonical  problem  -  that  of  a  line  source  over  an 
impedance  plane.  Then  the  theory  is  applied  to  the  surface  wave  excited  bv  a 
line  source,  on  a  wedge  vjith  variable  surface  impedance.  The  result  agrees 
precisely  with  the  asymptotic  form  of  the  exact  solution.  Another  application 
is  made  to  the  surface  wave  excited  on  a  cylinder  bv  a  line  source.  The  result 
also  agrees  with  the  exact  solution. 
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1.  Introduction 

A  surface  wave  is  a  wave  which  travels  along  a  surface  with  its  field 
confined  to  the  neighborhood  of  the  surface.  It  differs  from  other  waves  with 
this  property  in  that  it  can  exist  by  itself  on  a  flat  surface  in  a  homogeneous 
medium.  The  evanescent  waves  produced  by  total  reflection  are  also  confined  to 
the  neighborhood  of  a  surface,  but  they  cannot  exist  without  the  incident  field. 
On  the  other  hand,  the  diffracted  waves  associated  with  the  usual  surface  rays 
can  exist  by  themselves  on  ciirved  surfaces,  but  not  on  flat  ones.  A  surface  wave 
is  characterized  by  a  complex  propagation  constant  k,  =  ol)/c,  ,  where  <J0  is  the 
angular  frequency  of  the  field  and  c-,  is  the  complex  velocity.  Either  k-  or 
c^  determines  the  phase  velocity  and  the  attenuation  rate  of  the  wave  along 
the  surface.  We  shall  assume  that  k,  and  c,  are  known  and  investigate  the 
manner  in  which  such  waves  are  excited,  how  they  propagate  and  the  rate  at 
which  their  field  diminishes  with  distance  from  the  surface.  We  shall  do  this 
by  constructing  a  geometrical  theory,  analogous  to  geometrical  optics  and  to  the 
geometrical  theory  of  diffraction  [^  J  ,  which  will  permit  us  to  visualize  tne 
excitation  process  and  to  calculate  the  amplitude  of  the  surface  wave  in  rather 

general  cases. 

la  the  next  three  sections  we  formulate  our  theory.  In  the  fifth  and 

sixth  sections  we  determine  the  excitation  and  radiation  coefficients  which 

enter  the  theory.  This  is  done  for  scalar  fields  which  satisfy  the  reduced 

wave  equation  and  satisfy  an  impedance  condition  on  the  surface.  The  problem  of 
designing  the  source  to  maximize  the  surface  wave  amplitude  is  also  treated. 
Then  the  method  is  used  to  determine  the  surface  wave  which  a  line  source  excites 
on  a  wedge  with  variable  surface  impedance.  The  resxilt  is  found  to  coincide 
with  the  high  frequency  asymptotic  form  of  the  exact  solution  of  this  problem. 
Finally  the  surface  wave  on  a  circular  cylinder  is  determined  by  omt   method 
and  also  found  to  coincide  with  the  asymptotic  form  of  the  exact  solution. 
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2.   Complex  rays 

Rays  in  space  can  be  defined  in  terms  of  the  propagation  velocity 
c(x)  by  Fermat's  principle  or  by  the  differential  equations  which  follow 
from  it.  Complex  rays  can  be  defined  by  apolying  Fermat's  principle  to 
complex  curves,  or  bv  considering  complex  solutions  of  the  ray  equations. 
In  either  case  it  is  necessary  that  c(x)  be  defined  for  complex  values  of 
X  ,  which  is  the  case  if  c(x)  is  an  analytic  function  of  x.  We  can  also 
define  complex  surface  rays  by  applying  Fermat's  principle  to  complex 
curves  on  the  surface  3  or  by  considering:  complex  solutions  of  the  equations 
for  surface  rays.  In  this  case  it  is  necessary  that  both  the  velocity 
c^  (x)  and  the  surface  S  be  defined  for  complex  x  . 

Next,  by  means  of  Fermat's  principle,  we  may  define  a  ray  joining  two 
points  Q  and  P  in  space  and  having  an  arc  on  S.  It  follows  from  Fermat's 
principle  that  such  a  ray  consists  of  a  complex  space  ray  from  Q  to  some 
point  Q^  on  S,  a  complex  surface  ray  from  Q_  to  another  point  P  on  S,  and 
a  complex  space  ray  from  P.  to  P.  In  addition  it  follows  that  at  both  Q^ 
and  P  the  unit  tangents  I  and  T  to  the  space  and  surface  rays  respectively 
must  satisfy  the  condition 

(1)  c""'"I  +  c"-^T  »  (c'-^I«>J+c^-'-T«N)N. 

Here  N  denotes  the  unit  normal  to  S  at  Q,  or  P^.  The  vectors  are  oriented  such 
that  -I  and  T  both  point  in  the  direction  of  increasing  arclength  along  the  ray. 
See  Figure  1.  Equation  (l)  is  equivalent .to  the  geometrical  conditions  that  I 
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FIGURE    I 


The  complex  ray  from  Q  reaches  the  surface  at  the  complex  point  Q^   where 
it  produces  a  complex  surface  ray.     This  ray  travels  along  the  surface  with 
the   velocity  c,   and  sheds  complex  space  rays.     The  ray  which  leaves  the  surface 
at  P^   passes  through  P.     The  angles  6  and  ©,    are  determined  by  Snell's  law  (l). 


3  - 


and  T  be  coplanar  with  N  and  lie  on  opposite  sides  of  N,  and  that  the  angle 
©  between  I  and  N  satisfy  Snell's  law 


(2) 


sin  ©  =  — 


The  preceding  considerations  enable  us  to  determine  all  the  ray  paths  which 
contain  arcs  on  S  along  which  surface  wave  propagation  occurs. 

3»  Amplitude  and  phase 

The  field  on  any  ray  is  composed  of  a  phase  and  an  amplitude.  On  a 
surface  ray  from  Q^  to  P  the  complex  phase  change  0  ((i.,P. )  is  the  follow- 
ing integral  with  respect  to  arc  length  s 


(3) 


On  a  space  ray 


?^(Ql,Pl) 


\     °iP^^^]  "  { 


k^  x(s)  ds  . 


(U) 


p         p 

J         c|_x(s)J     > 


1  -   ^    1 

The  amplitude  variation  along  anv  ray  deoends  upon  the  family  of  rays  in 
which  the  ray  under  consideration  is  continued.  For  a  complex  ray  the 
amplitudes  A  (Q  )  and  A  (P), which  are  vectors  for  vector  fields,  are 
related  by 


(?) 


A^(P^)  -  A^(Q^) 


d(y(Q^)c3^(P^) 
daCP^;c^(Q^; 


1/2 


-  u  - 


In  (5)  dotQ, )/dff(P^)  denotes  the  ratio  of  the  width  at  Q.  to  that  at  P. 
of  an  infinitesimal  band  of  complex  surface  containing  the  ray  from  Q^  to 
P^.  The  analogue  of  (5)  for  the  space  ray  from  P  to  P  is 

1/2 


(6) 


A(P)  -  A(P^) 


dr(p^)c(p) 


dr  (P)c(p^; 


In  (6)  dT  (P  )/dT  (P)  is  the  ratio  of  the  cross  sectional  area  at  P^ 
to  that  at  P  of  an  infinitesimal  tube  of  complex  rays  containing  the 
ray  from  P,  to  P.   A  similar  relation  holds  on  any  space  ray. 

The  amplitude  on  the  surface  ray  at  Q.  is  proportional  to  that  on 
the  space  ray  at  Q^  so  we  have 


(7) 


A  (Q^)  =  E(Q^)A(Q^) 


The  proportionality  factor  E(CI  )  in  (7)  is  an  excitation  coefficient. 
In  the  case  of  vector  fields,  it  is  a  matrix.   Similarly,  the  amplitude 
on  the  space  rayatP^is  proportional  to  that  on  the  surface  ray  there 


so  we  have 


C^)  A(P^)  =  R(P^)A^(P^)  . 

The  proportionality  factor  or  matrix,  R(P-,  )  in  (8)  is  a  radiation  coefficient. 

Upon  combining  the  above  results,  we  find  that  the  field  u(P)  at  P  is  given 

by 

1/2 


(9)  u(P)  -  R(P^)E(Q^)A(Q^) 


d<y(Q^)c^(P^)dr(P^)c(P) 


d(y(P^)c^(Q^)dr  (pjc(Pt) 


ip(Q^)+0l(<\,Pl)*0(Pl,P)'| 


Equation  (9)  expresses  the  field  at  P  on  a  ray  which  has  traveled  along  an 
arc  of  S  as  a  comolex  surface  ray  associated  with  surface  wave  propagation. 
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It  involves  the  incident  field  A(Q- )e 


15^(0^) 


,  various  geometrical  quantities, 
and  the  excitation  and  radiation  coefficients. 

We  shall  now  restrict  our  considerations  to  the  two  dimensional  case. 
In  addition  we  shall  assume  that  c  is  constant  so  that  the  space  rays  are 
straight  lines.  Due  to  the  two  dimensionality,  the  ratio  d(j(Q- )/dcr(P,  )  -  1. 
Because  the  rays  are  straight  we  have  0(P,,P)  =  kd(P^,P)  where  d(P^, P)  is  the 
complex  distance  from  P^  to  P.  In  addition  we  find  that 


(10) 


___Li- -  ri+(X+9)d(P,,?)sec9l 


-1 

P, 


In  (10)  X(P-,  )  denotes  the  curvature  of  S  at  P^,  ©(P,  )  is  the  solution  of 
(2)  at  P^,  and  9  is  the  derivative  of  6  with  respect  to  arclength  along 
the  ray  on  S,  evaluated  at  P^ ,  With  these  simplifications,  (9)  becomes 


(ll)u(P)  -  R(P^)E(Q^)A(Q^) 


c,(P,) 


h^ 


1/2 


l+C^+e)d(P  ,P)sece 


-1/2 


p(Q^)+0^(Q^,P^)+kd(P^,P; 


h.       The  reciprocity  principle 

To  determine  the  coefficients  R(P, )  and  E(Q. )   we  will  now  apply  (11) 
to  the  case  in  which  the  field  is  produced  by  an  isotropic  line  source, 
represented  by  the  point  Q  in  two  dimensions.     Then  the  incident  field  is 
given  by 


(12) 


A(Q^) 


^^^^^       .    r.fn.  o  M-1/2  i''^^'^'^) 


A^[d(Q,Q^)j' 


In  (12)  A  is  a  constant  characteristic  of  the  source,  with  the  dimensions 
^     o  ' 

1/2 
of  (length)  '   times  those  of  u.  When  use  is  made  of  (12),  (11)  becomes 
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-1/2 


(13)     u(P)  =  R(P;L)E(q^)A^|^d(Q,Q^)j"-'/''  l+U+e)d(P^,P)sece 


-1/2 

P, 


c,(P,) 


^I^ 


1/2  , 
e 


irkd(Q,Q^)+0(Q^,P^)+kd(P^,P)1 


We  now  make  use  of  the  reciprocity  principle  which  requires  that  a 
source  at  P  must  preduce  at  Q  the  same  field  as  an  equal  source  at  Q 
would  produce  at  P,  This  implies  that  the  right  side  of  (13)  must  be 
unchanged  when  P  is  interchanged  with  Q  and  P^  with  Q, ,  Frcm  this  con- 
stancy we  obtain  the  condition 


(Hi)         R(P^)E(Q^)[d(Q,n^)-| 


-1/2 


l+(X+9M(P3^,P)sec9 


-1/2 
P, 


J  1 


c,(P,) 


\^ 


1/2 


R(Q^)E(P^)[d(P,P^)] 


-1/2 


l+(Xt©)d(Q^,Q)sec© 


,1-1/2  pl(«l)' 


1/2 


Eouation  (1)4)  enables  us  to  express  the  excitation  coefficient  E  in  terms  of 
the  radiation  coefficient  R  by  the  relation 


(15) 


;(Q^)  =  a^R(Q^)  U...(X+e)sece  ■•'■/^  c^(Q^) 


10^ 


Here  a  is  a  constant  with  dimensions  (time)X(length)  '  •  From  (15)  we  see  that 
E(Q.)  depends  nponj^."   l/d(Q,Q.),  the  curvature  of  the  incident  wavefront 
at  Q. ,  in  addition  to  properties  of  the  surface  at  Q. .  It  will  be  conven- 
ient to  introduce  R  (Q^)  =  a  c^ (Q- )    R(Q, )  in  order  to  avoid  carrying  along  the 
factor  a. 
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If  we  now  use  (13')  to  eliminate  E(Q^)  from  (13),  we  obtain 


(16)   u(P)  =  r'(P^)r'(Q^)A^  l+(^+5)d(Q^,Q)sec9  "^^^  l+(X+5)d(P^,P)sec9 


1 


-1/2 

P, 


X  e 


irkd(Q,Q^)+0(Q^,P^)+kd(P^,P)1 


This  expression  for  u(P)  applies  in  the  two  dimensional  case  when  c  is 
constant  and  the  field  is  produced  by  a  line  source  at  Q.  For  a  more 
general  incident  field  we  must  use  (1^)   in  (11). 

If  the  surface  S  is  a  plane  and  if  c^  is  constant  then  p  is  infinite, 

§  is  zero,  r'(P  )  =  r'(Q^)  =  R',  0(Q^,P^)  »  k^d(Q^,P^)  and  (16)  becomes 

,-   irkd(Q,Q.)+k.d(a,P  )*kd(PT,P)'I 
(17)        u(P)  '   R  ^A^e  L     1111      1   J    , 

^,       Determination  of  the  coefficients 

To  complete  the  prescription  for  calculating  the  field,  we  must 
determine  the  radiation  coefficient  R  .  Since  R  is  dimensionless,  R 
has  the  dimensions  of  (length)"  ''^''.  Therefore  it  will  be  mor-e  conven- 
lent  to  consider  the  dimensionless  product  k    R  (P,).  This  product 
must  depend  upon  the  type  of  field  under  consideration.  We  shall  assume 
that  it  depends  upon  ©(P^)  and  possibly  upon  other  physical  constants  of 
the  medium  and  the  surface  at  P^,  but  not  upon  any  geometrical  properties 
of  the  surface  at  P^ •  Therefore  it  can  be  found  by  considering  the  special 
problem  of  a  line  source  in  the  presence  of  a  uniform  plane  surface  in  a 
homogeneous  medium.  By  comparing  the  exact  solution  of  such  a  problem  with 
(17),  after  expanding  the  exact  solution  for  large  co  (or  large  k  and  k^), 
R  can  be  determined.  In  order  to  facilitate  such  a  comparison,  we  will 
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now  calculate  the  distances  which  appear  in  (1?)  in  terras  of  the  coordinates 
of  P  and  Q. 

Let  us  assume  that  the  surface  occupies  the  plane  y  »  0,  that  the  source 
Q  is  located  at  (0,y^),  and  that  the  point  P  is  at  (x,y)  with  x  >  0,  The 
ray  from  Q  to  Q^  is  a  straight  line  which  makes  with  the  v-axis  the  angle  9, 
determined  by  (2).  Similarly  the  ray  from  P.  to  P  is  a  straight  line  making 
the  same  angle  with  the  y-axis.  Therefore  the  coordinates  of  Q^  are 
(y^tanP,0)  and  those  of  P^  are  (x-3rtan9,0).  See  Figure  2.  Then  we  find  that 

'i(Q»Q]_)  =  y^secG,  d(P^P)  =  ysecO  and  d(Q^,P^)  =  x-(y^+y)tan9.  When  these 

distances  are  inserted  into  (17),  it  becomes 

,„   ik  lxlsin9+(y  +y)cos9 


(18)  u(x,y)  -  R 


0 


Here  x  has  been  replaced  by  |x|  in  order  that  this  result  remain  valid 
for  X  <  0.  In  obtaining  (18)  from  (17)  we  used  the  relation  k/k,  •»  sin9, 
which  follows  from  (2). 

In  the  next  section  we  shall  obtain  an  exact  solution  of  the  problem 
just  treated  for  a  scalar  field  which  satisfies  the  reduced  wave  equation 

(19)  (V^+k^)u  -  5(x)5(y-y^)  ,       y>0. 

On  the  plane  y  =  0,   u  will  satisfy  the  impedance  boundary  condition 

(20)  u     -  -ikZu  . 

In  (20)  the  complex  constant  2  is  the  surface  impedance.  In  order  for 
the  surface  to  support  a  surface  wave,  it  is  necessary  that  Z  satisfy 
the  conditions 
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H 


NO  SURFACE  WAVE 


SURFACE    WAVE 

P(x,y) 


/   ///  . 


/    /d^/     ///////     p     /////     /^ 


^'       (-yo/v5)mZ,o) 

(0,-yo) 


FIGURE  2 


The  surface  wave  produced  by  a  line  source  at  Q  above  a  plane  of  surface 
impedance  Z.  The  figure  shows  a  ray  from  Q  which  makes  the  complex  angle  © 
with  the  normal  to  the  plane.  The  angle  ©  is  determined  by  equation  (22). 
This  ray  gives  rise  to  a  complex  surface  ray  which  sheds  complex  space  rays. 
One  such  ray  is  shown  leaving  the  surface  at  P,  and  passing  through  P.  These 
rays  cover  the  region  to  the  right  of  the  real  line  (y  +  y  )x~"^  =-ImZ  provided 
Re  Z  =  0.  The  figure  is  symmetric  about  the  y  axis. 
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(21)  ReZ  >  0  ,  ImZ   <  0  . 

If  (18)  is  inserted  into  (20)  it  vields 

(22)  cos9  =-Z  . 

From  (22)  and  (2)  the  propagation  c^  is  determined  by  the  surface 
impedance  through  the  relation 

(23)  c^  »  c(l-Z^)"^/^  . 

When  the  outgoing  solution  of  (19)  and  (20)  is  expanded  asymptotically 
for  large  values  of  k,  it  contains  a  surface  wave  of  thefbrm  (18)  with 
9  given  by  (22).  This  surface  wave  becomes  identical  with  (18)  if  in 
(18)  we  set 

(2li)  A^  =  (8nk)"^/2  e"^"/^ 

(25)  r'  =  (Snk)^/^  e^"/^  (cot9)^/^  . 

This  equation  determines  the  radiation  coefficient. 

Let  us  now  insert  (2)i)  and  (25)  into  (16).  This  yields  the  final 
result  for  the  surface  wave  excited  on  an  impedance  surface  by  an 
isotropic  line  source  of  amplitude  A  in  a  homogeneous  medium.  It  is 


(26)  u(P)  =  |tan©+(>l4)d(Q^,Q)secetanep^rtane+(>{+§)d(P^,P)secQtan9l'-'-'^^ 


X  e 


i[kd(Q^,Q^)+!Zf(Q^,P^)+kd(P^,P)J 
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If  the  line  source  is  not  isotropic  we  may  write  its  amplitude  as 
A  f(0).  Then  the  right  side  of  (26)  must  be  multiplied  by  f(©).  Here 
f (©)  represents  the  radiation  pattern  of  the  source.  It  is  proportional 
to  the  amplitude  on  a  ray  which  makes  the  angle  9  with  the  normal  to  the 
surface  3.  Usually  the  pattern  is  expressed  as  a  function  F(T)  of  the  angle 
7  which  the  ray  makes  with  a  fixed  direction.  By  determining  T  as  a 
function  of  9  from  the  shape  of   3  and  the  location  of  Q,  we  then  have 
f(9)  =  f[t(9)].  To  maximize  the  amplitude  of  the  surface  wave  it  is  nec- 
essary to  maximize  the  pattern  f(9)  in  the  complex  direction  9  of  the 
ray  which  excites  the  sxirface  wave.  This  directio,n  is  determined  by  (2) 
or  (22).  If  S  is  a  plane  and  if  Y  is  measured  from  the  line  through  Q 
normal  to  3,  then  T  -  9  and  therefore  f(9)  =  F(9). 

If  the  incident  field  is  arbitrary,  the  surface  wave  is  given  by 
(11),  (15)  and  (2^)*     Upon  combining  these  equations  we  obtain 


(27)   u(P)  -  (8nk)^/^u.(Q,)U'.tan9+(/<4'>gec9tan9|  ■'•/^rtan9t(>i:+9)d(P  ,P)sec9tan9| 

i[n/2+0^(QpP^)+kd(P^,P)J 


X  e 


Here  u. (Q)   =  A(Q)e  denotes  the  incident  field  at  Q^, 

6.       Siirface  wave  on  an  impedance  plane 

We  shall  now  solve  the  boundary  value  problem   (19)  and   (20),    also 
requiring  u  to  satisfy  the  radiation  condition.     The  solution  is 

■i      (T')  ^      f  r.i   ni  r,       ikrxcos3f+(y+y  )sin01 

(28)        u(x,y)=-^jj(l'to).5i    j|i^     e      L         '''='0         ''J,^. 
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In  (28)  r  denotes  the  distance  from  (0,y  )  and  C  is  a  contour  in  the  complex 
plane.   It  runs  from  n  -  ioo  to  n   along  the  line  Re0  »  n,  then  from  n  to  0 
along  the  line  Im0  =  0  and-  finally  from  0  to  ioD  along  Re^  =  0.  See  figure  3- 

To  obtain  the  asymptotic  form  of  u  for  large  k,  we  can  first  use 
the  asymototic  form  of  H^  '^(kr)  for  large  values  of  the  argument.  This 
yields  the  asymptotic  form  of  the  incident  field.  The  asymptotic  form 
of  the  reflected  field  is  obtained  by  performing  a  saddle  point  evaluation 
of  the  integral  in  (28).  The  saddle  point  0  which  lies  in  the  interval 
0  <  0  <  n  is  given  by 

(29)  0„  =  tan  ^ ^ 

|x| 

The  appropriate  contour  through  the  saddle  point  is  obtained  by  setting 
the  imaginary  part  of  the  phase  xcosi^+(y+y  )s±n0   equal  to  its  value  at 
^  ~  ^L*  I^  "s  Is't  ^  "   a+i3>  we  find  in  this  way  two  possible  contours 

(30)  a  «  5f  +  cos  (sechp)  . 

In  order  that  the  exponential  in  the  integrand  of  (28)  vanish  at  infinity 
on  the  contour,  we  must  choose 

(31)  a  =  0  +co3~  (sechB)        0  <  0 

(32)  a  =  0  -cos  (sechR)        3  >  0 

With  this  choice  of  signs,  the  contour  C  can  be  deformed  into  the  saddle 

contoTjT,  See  Figure  3.  Then  the  saddle  point  evaluation  of  the  integral 
yields  the  asymptotic  form  of  the  reflected  field. 


-  Ua  - 
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FIGURE   3 


Saddle  point  contoiir  used  in  the  asymptotic  evaluation  of  the 
reflected  field  for  a  line  source  above  a  plane  with  a  constant 
surface  impedance. 
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In  the  process  of  deforming  C  into  the  saddle  contour  it  may  cross 
the  Dole  of  the  integrand  which  is  at  0  given  by 


(33)  sin0^  =  -Z 


If  this  occurs,  the  residue  of  the  integral  at  this  pole  must  be  added 
to  the  saddle  point  contribution  to  u.  The  residue  is  just  the  surface 
wave.  It  is  easil^^  found  to  be  given  by 


2   ikyxl  /l^^  -  (y+yQ)z] 


(3U) 


We  have  already  stated  that  this  coincides  with  (18)  when  A  and  R  are 

given  by  (2I4)  and  (2^),  and  9  by  (22). 

To  determine  when  the  surface  wave  (3h)  is  present,  let  us  consider 
the  saddle  contoiir  which  passes  through  the  pole,  ^et  f     denote  the  value 
of  0  for  this  contour.  Because  Z  satisfies  (21),  the  pole  lies  in  the  upper 
half  plane  p  >  0.  Thus  (32)  yields  for  ^  the  resxilt 

(35)  ^  '   c.+cos~  (sechB-,)   . 

If  0  <  ?r ,  the  pole  is  crossed  and  the  surface  wave  occurs,  if  9^-  >  ^"^ 
it  does  not  occur.  From  (29)  we  thus  find  that  the  surface  wave  occurs 
if  and  only  if  the  following  condition  is  satisfied 

^3^^ <  tanra,+cos"-'-(s6ch0T)|   . 

In  (36)  a.   and  p.  are  the  real  and  imagipgry  parts  of  0.,  the  solution 
of  (33).  In  the  special  case  in  which  ReZ  -  0,  (33)  shows  that  a^  «  0. 
Then  p  =  -sinh"  (ImZ)  and  (36)  becomes 

(37)  IZ2.  <  •  ImZ  . 

See  figure  2. 
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7»   Surface  wave  on  a  wedge  with  a  non-uniform  impedance 

Let  us  now  apply  our  method  to  determine  the  surface  wave  which  a 
line  source  excites  on  a  wedge  with  a  non-uniform  surface  impedance.  If 
the  line  source  is  parallel  to  the  edge  of  the  wedge, this  problem  is 
two  dimensional.  Then  the  point  Q  at  which  the  source  is  located  can  be 
described  by  its  polar  coordinates  (p',?^  )  and  any  other  point  P  by  (o,0). 
The  origin  is  at  the  tip  of  the  wedge,  the  surfaces  of  which  are  at  j^  =  0 
and  ^  ■  a.  See  Figure  U.  If  the  irrci'^'^nt  wave  has  the  amplitude  (2[|)  then 
the  surface  wave  which  it  produces  on  either  surface  of  the  wedge  is  given  by  (26). 
To  determine  this  wave  excited  on  0  =  0  we  assume  that  Z,  the  surface 
impedance  of  the  wedge  on  0'  =  0,  is  proportional  to  p   .  Thus 

(38^  Z  -  Z^p-^ 

Here  Z  is  a  constant  which  satisfies  (21).  Then  from  (23) 
o 

(39)  c^  "  pc(p^-Z^)-^/2  . 

To  write  (26)  more  explicitly  we  must  evaluate  the  various  geometrical 
quantities  which  aopear  in  it.  We  .''hall  first  give  them  when  P  is  to  the 
right  of  Q,  so  that  pcos0  >  p'cos0  ,  Then 

(liO)  d(Q,Q^)  -  p'sin0  sec© 

(m)  d(P,P-)  -  psin^sec© 

(U2)         p(Qn)  =  p'cos0  +p»sin0  tan©  -  p'sec©  cos(9('  -©  ) 


-  13a  - 


0{P'4>) 


/  /  / / / // //////  /  /  /////// 


11 


-►  X 


J, 


<^  =  0 


* 


FIGURE    4 


Excitation  of  a  sixrface  wave  on  a  wedge  by  a  line  source  at  Q.  The 
coordinates  of  Q  are  (p',^).  A  complex  ray  from  Q  excites  a  surface  ray  at 
the  complex  point  Q-..  This  ray  travels  along  the  surface  shedding  complex 
space  rays.  The  ray  shed  at  P^  passes  through  P  which  has  coordinates  (p,j20. 
The  wedge  faces  are  at  ^  »  0  and  0  =  a.  The  complex  angles  ©  and  6,  are 
determined  by  equation  (22)  in  which  Z  is  given  by  equation  (38). 
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(l43)  d(P,  )  *  pcos0  -  p3in0tane  =  psec©cos(0+9) . 

From  (22)  and  (33)  we  have 

(U)4)  p(Q^)  =  -Z^sece' 

V6)  p(P^)  =  -Z^secP  . 

Now  from  (I42)  and  ()4)i)  we  determine  9  ,  and  from  (1|3)  and  ('45)  we  find  9, 

(h6)  e'  =  !Zf'+cc3"-'-(-Z^/p') 

(hi)  9  =  -5f^-co3'^(-Z^/p)  . 

i^rom  (3)   and   (39)   we  obtain 

P(P-,) 


(U8) 


p(Ql) 


k  pcos(^+9)tan©-p»cos(;2f  -9   )tan9   +Z    (9  -9) 


Upon  combining    ('4O),    ()4l)   and   ('48)   and  simplifying,   we  obtain  for  the   phase 
in   (26) 

('i9)  kd(Q,Q^)+?f^(Qj^,P^)+kd(P^,P) 

-   k[p(l-z2p-2)l/2  _p,(i.z2p,-2)l/2_2^^^^^') 

-Z   cos"-'-iZ    /p)+Z  cos"-^(Z    /p«)        . 
0  00  °         J 

To  calculate  the  amplitude  in  (26)  we  determine  9  and  5  from  (22) 
and  (38)  and  observe  that  >C  =  0.  First  we  find 

(50)  9'sec9'tan9'  »  9sec9tan9  -  -Z'^ 
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(?1)  tan9t-esecetaned(P,P^)  »  tan  co3'-'-(-Z^p"-'-)-(^  -  Z^-^-psinjZfsecfcos'-^C-Z  p""^)-^ 


.(p^Z^^-l)^/2 


(52) 


tan9*+e'3ece'tane'd(Q,Q^)  =  -(p'^z'^-l)-^^^ 


When  (5l),  (52)  and  ([i9)  are  used  in  (26),  the  final  result  for  u  when 
pcos^  >  p^co3^  becomes 

(53)     u(p,!^)  =  (p2z;2-i)-lA(p,2,-2.2)-lA 

ik[p(l-z2p-2)l/2.p,(i.z2p,-2)l/2_2^(j2j^5^') 

-Z  cos"''-(Z  /p)+Z  cos"^(Z  /p') 


><  exp  . 


>  • 


If  pcos/  <  p'cosf^  ,  u(p,9f)  is  obtained  by  interchanging  primed  and  unprimed 

quantities  in  (53).   In  this  case  if  P  is  near  the  tip,  then  p  <  <  |Z  |  and 

—1  r  —1    r  ?  ?i  1 /?  "1 

(53)  can  be  simplified.  First  we  replace  cos  (Z  /p)  by  ilog  p  (2  +|Z  -p  J   )  . 

Then  upon  omitting  some  small  terms,  (53)  becomes 


(5h) 


u(p,^)  =  e--A(p.2z-2.i)-^A 


exp, 


ik[-iZ^-p'(l-z2p.-2)V2.2^(^,j2f') 
+Z^cos"^(zyp«  )+iZ^log(2Zyp)J 


If  Q  is  also  near  the  tip  (5'i)  can  be  further  simplified. 

Let  us  now  compare  our  results  (53)  and  (5'j)  with  the  asymptotic  form 
of  the  exact  solution  of  the  corresoonding  boundary  value  problem.  We  seek 
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a  solution  u  of  (19)  with  the  source  at  (p',?f  ).  This  solution  must  be 
defined  in  the  region  outside  the  wedge^i.e.^in  the  region  0  <  (^  <  a.  On 
the  surface  ^  =  0  it  must  satisfy  (20)  with  Z  given  by  (38).  On  5^  =  a  we 
require  it  to  satisfy  u^v  =  0.  The  solution  must  also  satisfy  the  radiation 
condition. 

r2i 

The  solution  of  this  problem  has  been  found  by  L.  Felsen  -  -■  .  It 
contains  a  surface  wave  given  by 

(55)  u(o,!Z()  =  -in  coshn(a-gf)cosh.(a-g(')  ^    (kp  )H^J^  (kp  )   . 

a-i(kZ^)'^sinh^/ia      "^"^    <  "^"^    "" 

Here  p  and  p  denote  respectively  the  smaller  and  larger  of  the  two  quantities 
p'  and  p.  The  constant  -rj  is  a  solution  of  the  equation 

(56)  coth  a  =  -?i(kZ  )"""•    . 

For  large  k  the  solution  of  (56)  becomes 

(57)  n~ikZ^   • 

Thus    (55)  becomes 

ikZ   (i2f+?f')  (^. 

(58)  u(p,5f)  ~  nkZ^  e       °  '^kZ   ^^^<^"kZ^'"'>"'    • 

o  o 

We  now  insert  the  Debye  asymptotic  forms  of  the  Bessel  and  Hankel 
functions  for  large  k  with  order  less  than  the  argument  in  absolute  value. 
Then  (58)  becomes,  if  p  >  p'  >  |z  | 


(?■»    u(.,0)  ~  Z(,'Ui)-'-\/.zl)-^/''   .-iX^tfV) 


O 

X  exp. 


ik[(p'2-z2)^/2-(p2-z2)V2_2^^^^-l(2yp,)^2^^^^-l(2yp)j 
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This  agrees  exactly  vrith  (53).   If  p  <-  p'  we  must  interchange  p  and  p'  in 

(59)  and  the  result  again  agrees  with  (53)  when  it  is  similarly  changed. 
If  D  <  <  !Z  I  and  d'  >  |Z  I  then  (56)  becomes  identical  with  (5U)  when  we 
insert  the  asymptotic  form  of  the  Bessel  function  for  large  k  with  order 
larger  in  magnitude  than  the  argument,  and  omit  some  small  terns.  The 
agreement  between  the  asymptotic  form  of  the  exact  solution  {'^S)   and  the 
solution  (53)  or  (510  given  by  our  'aethod,  is  a  confirmation  of  our  theory. 

8.   Surface  wave  on  a  cylinder 

Let  us  finally  use  our  method  to  find  the  surface  wave  excited  on  a 
circular  cylinder  of  radius  a  by  a  line  source  parallel  to  the  cylinder 
axis.   If  the  amplitude  of  the  wave  from  the  source  is  given  by  (2li)  then 
the  surface  wave  is  given  by  (26).  We  assume  that  the  cylinder  has  a  constant 
surface  impedance.   Then  Q  in  (26)  is  a  constant  determined  by  (22).  Let  the 
polsu:  coordinates  of  Q  and  P  be  (p',0)  and  (p,^  respectively.  There  are  two 
rays  from  Q  and  P  which  traverse  the  cylinder  in  opposite  directions.  We  first 
consider  the  shorter  ray.  See  Figure  5.  The  two  corresponding  pairs  of  points 
Q- ,  P^  are  foimd  to  have  the  coordinates  p(Q-,)  =  p(P-,)  ■  a  and 

(60)  5^(0^)  »  cos'^(-Z)-sin"-'-(|,  V^^  )  . 


(61) 


Then  we  find  that 


^(P  )  =  0-cos'^(-Z)+sin"^(-  S^   ) 


(62)  d(0,a,^)  '   p' 

(63)  d(P,P^)  =  p 


^/-    .2.  .  ,a 


^.)'(1-Z^).(^-,)Z 


1 

^)Z 

p       p  J 


/-(i)  (i-z^)  .  (L) 
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FIGURE    5 


Excitation  of  a  surface  wave  on  a  circular  cylinder  by  a  line  source  at  Q. 
Two  complex  rays  from  Q  hit  the  cylinder  at  the  complex  points  Q,  and  Qp  where 
they  make  the  complex  angle  ©  with  the  normal.  This  angle  is  determined  by 
equation  (22).  The  two  complex  surface  rays  which  they  excite  travel  around 
the  cylinder  in  opposite  directions.  The  complex  space  rays  which  they  shed 
at  P-  and  Pp  pass  through  P. 
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(6l4)     0^{Q^,?^)  '  k^ap(P^)-(2f(Q^)-| 

=  ka  j/l^  r0-2cos"^(-Z)+sin"^f-v/i^  )+sin-ri    /^)1    . 

Similar  results  hold  for  the  longer  ray  except  that   ^  is  replaced  by  2tt-0 

in   (61)  and   (6)4).     Upon  inserting  these  expressions  into    (26)   a-;d  simplifying, 

we  obtain 


(65)    viioj) 


l^[(P/-(l.z2)]"'/^[;£')2.a.z2)] 


/:? 


-lA 


X       63 


xp  ika    e  /l-(^^(1-Z^)   ^  /l-(f   )^(l-zF) 


a  •       p 


P» 


+  vA^Jsin"^fi    v^C?')+sin"-(^       /T?) 
I  P  P' 

+  2Z-2     v/l-7  cos"-'"(-Z) 


X     exp  ika    \/i-Z     0+exp  ika    \/i-Z      (2n-0)     • 

There  are  also  other  rays  from  Q  to  P  each  of  which  encircles  the  cylinder 

n  times  before  leaving  at  P..     The  phase  0^(Q^,P^)   for  such  rays  exceeds 

2  1/2 
that  given  by  (6U)  by  the  amount  2nnk-a  =  2nnka(l-Z    )        .Their  fields  are  com- 

Duted  from   (65)   by  including  this  additional  phase.     When  the  fields   on 

these  rays  are  added  together  at  P,    a  geometric  series  results.      Its  sum 

is  the  total  surface  wave  at  P  which  is 
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(66)   „(p,^)=  ^[(£)^.(i.z2)]-iA[(j')2-(i-z2)]-iA 


x   exp  ika 


§  /i-(^)''(i-z==^)  ^  I'  yC(i,)'(i.z^) 


V^L^Jsin"^(-  /^)+sin"^(|,  v^)l 


+  2Z  -  2  /l-Z"^  cos"-'-(-Z) 


exp 


ika  v^-Z^0+exp  ika  /l-Z^(2n-0)  i 


l-expUnika  v^l-?J 

We  shall  now  compare  this  result  with  the  solution  of  the  corresponding 
boundary  value  problem.  That  problem  was  formulated  and  solved  in  [3]  for 
a  cylinder  with  an  impedance  Z.  However,  Z  was  not  assumed  to  satisfy  (21) 

so  that  surface  waves  of  the  type  we  are  considering  were  not  found.  Instead 
only  the  diffracted  fields  associated  with  the  real  surface  rays  were  con- 
sidered. Those  fields  were  obtained  by  finding  the  solutions  v  of  the  equation* 


(67) 


H^^^  (ka)  -  -iZH^^^ka) 


Each  root  v  yielded  a  wave  given  by  [see  reference  [3],  p  U3,  eq.(Li)]| 


,(2) 


(68 


^„  r—  cosv  (e-n)    oHv_  (ka)    ,    .  ,.. 

^    t         ^i^^m"    l^nHi^^(ka)   ^'" 


where  i\  F(x)  =■  ^  F(x)+iZF(x), 
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When  Z  satisfies  (21)  the  equation  (6?)  has  an  additional  solution 
which  corresponds  to  the  surface  wave.  It  can  be  found  by  using  in  (6?) 
the  Debye  asymptotic  fonnula  for  the  Hankel  function  for  order  and  argu- 
ment both  large  and  |v|  >  ka.  Then  (6?)  becomes 


(69)  i  /I 


[v2-(ka)2] 


1/2 


"  [v2-(ka)2]lA 


-  iZ 


}  exp 


r  2  ,,   s2nl/2      ^-1  V 

[v  -(ka)  I  •  +vcosh  r— 


\,   <=  0 


The  solution  of  (69)  is 


(70) 


V  =  ka(l-Z^)-'-/^ 


IVhen   (70)   is  used  in   (68)   and  the  asymptotic  form  of  the   Hankel  function 
is  employed,    (68)   becomes 


(71)     u(p,{2f)   = 


A^ 


(|) -d-z'^) 


lA 


(J')^-(i-z^) 


lA 


X   exp  ika 


£     /l.(i)2(i.z2).Pl    v/l-(i,)2(i.z^).     /I? 


sin"^(i    yC?) 
P 


+sin 


■\t,  ^) 


,  +  2Z-2 


vC? 


i  -  icosh'-"-    /.-Z^ 


exp 


ika    -/l-Z     0+exp  ika     J\-Z     (2n-5^) 


l-exp[2nika  /l-Z*^  J 


This  is  the  asymptotic  form  of  the  exact  surface  wave  (68) .  After  some  algebraic 
manipxilations,  it  agrees  exactly  with  the  result  (66)  given  by  our  method. 
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Comnanding  Officer 
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Attn:  E.T.  Allen,  Library-  Supervisor 

Bjorksten  Research  Labs,  Inc. 

P.  0.  Box  265 

Madison,  Wisconsin 

Attn:  Mrs.  Fern  B.  Korsgard 


Convalr,  A  Division  of  General  Dynamics 

Corp. 
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Convair,  A  Division  of  General  Dynamics 
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151?  Industrial  Way 
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Dome  and  Margolin,  Inc. 
29  New  York  Avenue 
Kestbury,  Long  Island,  N.  I. 
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P.O.  Box  200 

Long  Beach  1,  California 

Attn:  Engineering  Library  (C-2S0) 

Douglas  Aircraft  Co.,  Inc. 
"27  Lapham  Street 
El  Segundo,  California 
Attn:  Engineering  Library 

Douglas  Aircraft  Company,  Inc. 
3000  Ocean  Park  Boulevard 
Santa  Monica,  California 
Attn:  P.T.  Cllne 

Eq.  Sec.  Reference  Files, 

Eq.  Eng.  A250 

Douglas  Aircraft  Companv,    Inc. 
jrino  North  Memorial  Drive 
Tulsa,    Oklahoma 
Attn:     Engineering  Library,   D-250 

Electronics  Communication,    Inc. 
1R30  York  Road 
Tiwonium,   Maryland 

Emerson  and  Cuming,    Inc. 
R69  Washington  Street 
Canton,   Mass. 
Attn:     Mr.  W.   Cuming 

Emerson  Electric  Kfg.   Co. 

8100  West  Florissant  Avenue 

St.   Louis  21,    Missouri 

Attn:     Mr.   E.R.   Breslln,    Librarian 

Sylvanla  Elec.    Prod.    Inc. 
Electronic  Defense  Laboratory 
P.O.  Box  205  -    'Unci) 
Mountain  View,    California 
Attn:     Library 

Falrchild  Aircraft  Division 
Fairchild  Eng.    and  Airplane  Corp. 
Hagerstown,   Maryland 
Attn:     Library 

Farnsworth  Electronics  Company 
3700  East  Pont  lac  Street 
Fort  Wayne  1,    Indiana 
Attn:     Technical  Library 

Federal  Telecommunication  Labs. 
500  Washington  Avenue 
Nutley  10,    New  Jersey 
Attn:     Technical  Library 

The  Gabriel  Electronics 
Division  of  the  Gabriel  Company 
13?  Crescent  Road 
Needham  Heights   9li,   Mass. 
Attn:     Mr.   Steven  Galagan 


General  Electric  Advanced  Electronics  Center 
Cornell  University 
Ithaca,   New  York 
Attn:      J.   B.   Travis 

General  Electric  Company 

Electronics  Park 

Syracuse,   New  York 

Attn:     Documents  Library,  B.  Fletcher 

Building  3-lli3A 

General   Precision  Laboratory,    Inc. 

63  Bedford  Road 

Pleasantvllle,    New  York 

Attn:     Mrs.  Mary  G.   Herbst,    Librarian 

Goodvear  Aircraft  Corp. 

1210  Masslllon  Road 

Akron  15,  Ohio 

Attn:  Library  D/120  Plant  A 

Granger  Associates 

Electronic  Systems 

966  Commercial  Street 

Palo  Alto,  California 

Attni  John  V.  N.  Granger,  President 

Grumman  Aircraft  Engineering  Corporation 
Bethpage,  Long  Island,  N.  Y. 
Attn:  Mrs.  A.  M.  Gray,  Librarian 

Engineering  Library,  Plant  No.  5 

The  Hallicrafters  Company 

UljOl  West  ?th  Avenue 

Chicago  2U,  Illinois 

Attn:  La Verne  LaGloia,  Librarian 

Hoffman  Laboratories,  Inc. 
3761  South  Hill  Street 
Los  Angeles  7,  California 
Attn:  Engineering  Library 

Hughes  Aircraft  Company 
Antenna  Department 
Microwave  Laboratory 
Building  12,  Room  2617 
Culver  City,  California 
Attn:  M.  D.  Adcock 

Hughes  Aircraft  Company 
Florence  and  Teale  Streets 
Culver  City,  California 

Attn:  Dr.  L.C.  Van  Atta,  Associate  Director 
Research  Labs. 

Hycon  Eastern,  Inc. 
75  Cambridge  Parkway 
Cambridge,  Mass. 
Attn:  Mrs.  Lois  Seulowltz 
Technical  Librarian 

International  Business  Machines  Corp. 

Military  Products  Division 

590  Madison  Avenue 

New  York  33,  New  York 

Attn:  Mr.  C.F.  McElwain,  General  Manager 

International  Business  Machines  Corp. 
Military  Products  Division 
Owego,  Mew  York 

Attn:  Mr.  D.  I.  Marr,  Librarian 
Department  U59 

International  Resistance  Company 
IjOl  N.  Broad  Street 
Philadelphia  8,  Pa, 
Attn:  Research  Library 

Jansky  and  Bailey,  Inc. 
1339  Wisconsin  Avenue,  N.  W. 
Washington  7,  D.  C. 
Attn:  Mj.  Delmer  C.  Ports 

Dr.  Henry  Jasik,  Consulting  Engineer 
298  Shames  Drive 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 

Electromagnetic  Research  Corporation 
711  llith  Street,  N.  W. 
Washington  5,  D.  C. 
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Lockheed  Aircraft  Corporation 

2555  N.  Hollywood  Way 

California  Division  Engineering  Library 

Department  72-75,  Plant  A-1,  Bldg.  63-1 

Burbank,  California 

Attn!   N.  C.  Harnois 

The  Martin  Companv 

P.  0.  Box  179 

Denver  1,  Colorado 

Attn:  Mr»  Jack  McCormick 

The  Glenn  L.  Martin  Company 
Baltimore  3,  Maryland 
Attn:  Engineering  Library 
Antenna  Design  Group 

Maryland  Electronic  Manufacturing  Corp. 
5009  Calvert  Road 
College  Park,  Maryland 
Attn:  Mr.  H.  Warren  Cooper 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,  Rhode  Island 

The  W.  L.  Maxson  Corporation 

b60  West  l!jth  Street 

New  York,  N.  T. 

Attn:  Miss  Dorothy  Clark 

McDonnell  Aircraft  Corporation 

Lambert  Saint-Louis  Municipal  Airport 

Box  516,  St.  Louis  3,  Missouri 

Attn:  R.  D,  Detrlch,  Engineering  Library 

McMillan  Laboratory,  Inc. 

Brownville  Avenue 

Ipswich,  Massachusetts 

Attn:  Security  Officer.  Document  Room 

Melpar,    Inc.  (2) 

3000  Arlington  Boulevard 

Falls  Church,    Virginia 

Attn:     Engineering  Technical  Library 

Microwave  Development  Laboratory 
90  Broad  Street 
Babson  Park  57,   Massachusetts 
Attn:     N.   Tucker,    General  Manager 

Microwave  RaiJiation  Comr'any  Inc. 

19223  South  Hamilton  Street 

Garden*,    California 

Attn:     Mr.   Morris   J.  Ehrlich,    Prpsident 

Chance  Vought  Aircraft,    Inc. 
931b  West   Jefferson  Street 
Dallas,   Texas 
Attn:     Mr.   H.   S.   White,    Librarian 

Northrop  Aircraft,    Inc. 
Hawthorne,    California 

Attn:     Mr.  E.   A.   Freltas,    Library  DeptJll45 
1001  E.  Broadway 

Remington  Rand  Univ.   -  Division  of  Sperrv 

Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,    Pennsylvania 
Attn:     Mr.    John  F.   McCarthy 

R  and  D  Sales  and  Contracts 

North  American  Aviation,    Inc. 
1221h  Lakewood  Boulevard 
Downey,   California 
Attn:     Engineering  Library  U95-115 

North  American  Aviation,    Inc. 
Los  Angeles   Intemati'^nal   Airport 
Los  Angeles  Ii5,    California 
Attn:     Engineering  Technical  File 

Page  Coimminicstlons  Engineers,    Inc. 
710  Foiirteenth  Street,    Northwest 
Washington  5,    D.   C. 
Attn:     Librarian 

Phllco  Corporation  Research  Division 

Branch  Library 

U700  Wlssachlckon  Avenue 

Philadelphia  UU,    Pa. 

Attn:     Mrs.  Dorothy  S.  Collins 


Pickard  and  Burns,    Inc. 

210  Highland  Avenue 

Needham  9lj,   Mass. 

Attn:     Dr.   J.  T.   DeBettencourt 

Polytechnic  Research  and  Development 

Company,    Inc. 

202  Tiilary  Street 

Brooklyn  1,    New  York 

Attn:     Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,    Mass. 
Attn:     Dr.    John  Ruze 

Radiation,    Inc. 

P.  n.   Drawer  37 

Melbourne,    Florida 

Attn:     Technical  Library,   Mr.  M.L.  Cox 

Radio  Corp.   of  America 

RCA  Laboratories 

Rocky  Point,    New  York 

Attn:      P.   S.   Carter,    Lab.  Library 

RCA  Laboratories 
David  Sarnoff  Research  Center 
Princeton,    New  Jersey 
Attn;      Miss  Fern  Cloak,   Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electronic  Proincts 
Building  10,   Floor  7 
Camden  2,   New  Jersey 
Attn:     Mr,   Harold   J.   Schrader 

Staff  Engineer,    Organization 

of  Chief  Technical 

Administrator 

The  Ramo-Wooldridge  Corporation 
P.O.   Box  Ij5'i53  Alrnort  Station 
Los  Angeles  Ii5,    California 
Attn:      Margaret  C.  Whitnah, 
Chief  Librarian 

Hoover  MiC'^owave  Co. 
9592  Baltimore  Avenue 
College  Park,    Marvland 

Director,    USAF  Project  RAND 

Via:     Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  "-'ain  Street 

Santa  Monica,  California 

Rantec  Corporation 

Calabasas,  California 

Attn:  Grace  Keener,  Office  Manager 

Raytheon  Manufacturing  Company 
Missile  Systems  Division 
Bedford,  Mass. 
Attn:  Mr.  Irving  Goldstein 

Raytheon  Manufacturing  Company 
Wayland  Laboratory,  State  Road 
Wavland,  Mass. 
Attn:  Mr.  Robert  Borts 

Raytheon  Manufacturing  Comreny 
Wavland  Laboratory 
Wayland,  Mass. 

Attn:  Miss  Alice  G.  Anderson, 
Librarian 

Republic  Aviation  Corporation 
Farmingdale,  Long  Island,  N.  Y. 
Attn:  Engineering  Library 

Thru:  Air  Force  Plant  Representative 
Republic  Aviation  Corp. 
Farmingdale,  Long  Island,  N.Y. 

Rheem  Manu-''acturlng  Coapany 
9236  East  Hall  Road 
Downey,  California 
Attn:  J.  C.  Joerger 

Trans-Tech,  Inc. 
P.  0.  Box  3U6 
Frederick,  Maryland 


Ryan  Aeronautical  Companv 
Lindbergh  Field 
San  Diego  12,  California 
Attn:  Librar/  -  unclassified 

Sage  Laboratories 
159  Linden  Street 
Wellesley  81,  Mass. 

.Sanders  Associates 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.  R.  Wild,  Library 

Sandia  Corporation,  Sandla  Base 

P.O.  Box  5''00,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroscope  Company 

Great  Neck,  Long  Island,  New  York 

Attn:  Florence  W.  Tumbull,  Engr.  Librarian 

Stanford  Research  Institute 

Menlo  Park,  California 

Attn:  Library,  Engineering  Division 

Sylvania  Electric  Products,  Inc. 
100  first  Avenue 
Waltham  5li,  Mass. 

Attn:  Charles  A.  Thornhill,  Report  Librarian 
Walthan  Laboratories  Library 

Systems  Laboratories  Corporation 
lli''52  Ventura  Boulevard 
Sherman  Oaks,  California 
Attn:  Donald  L.  Margerum 

TRG,  Inc. 

17  Union  Square  West 

New  York  3,  N.  Y. 

Attn:  M.  L.  Henderson,  Librarian 

A.  S.  Thomas,  Inc. 

161  Devonshire  Street 

Boston  10,  Mass. 

Attn:   A.  S.  Thomas,  President 

Bell  Telephone  Laboratories 
Murray  Hill 
New  Jersey 

Chu  Associates 
P.  0.  Box  3"? 
Whitcomb  Avenue 
Littleton,  Mass. 

Microwave  Associates,  Inc. 
Burlington,  Mass. 

Raytheon  Manufacturing  Companv 
Missile  Division 
Harlwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  Laboratory 
225  Crescent  Street 
Waltham,  Mass. 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  50li,  Sunnyvale,  California 
Attn:  Miss  Eva  Lou  Robertson, 
Chief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  W.  C.  Hoffman 

Commander 

AF  Office  of  Scientific  Research 

Air  Research  and  Development  Command 

llith  Street  and  Constitution  Avenue 

Washington,  D.  C. 

Attn:     Mr.  Ottlng,   SRY 

Westinghouse  Electric  Corp. 
Electronics  Division 
Friendship  Int'l  Airport  Box  7L6 
Baltimore  3,   Maryland 
Attn:     Engineering  Library 
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WheelRr  Laboratories,  Inc. 
122  Cu-cter  Hill  Road 
Great  Neck,  New  York 
Attn:  Mr.  Harold  A.  Wheeler 

Zenith  Plaatics  Co. 
Box  91 

Oardena,  California 
Attn;  Mr.  S.  S.  Oleesky 

Library  Geophysical  Institute 
of  the  University  of  Alaska 
College 
Alaska 

University  of  California 

Berkeley  U,   California 

Attn:  Dr.  Samuel  Silver, 

Prof.  Engineering  Science 
Division  of  Elec.  Eng.  Electronics 
Research  Lab. 

University  of  California 
Electronics  Research  Lab. 
332  Corv  Hall 
Berkeley  li,  California 
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University  of  Pennsylvania 

Institute  of  Cooperative  Research 

3UO0  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Dept.  of  Electrical  Engineering 

University  of  Tennessee 
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ATTN:  Dr.  Pullen  Keats 
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b935  Whitehaven  Way 
San  Diego,  California 

University  of  Minnesota 
Institute  of  Technology 
Minneapolis,  Minnesota 
Attn:  Dean  Athelston  Spilhaus 

Ohio  State  University 
Columbus,  Ohio 
Attn:  Prof.  C.  T.  Tal 
Department  of  Electrical  Eng* 

Naval  Research  Laboratories 

Washington  25,  D.  C. 
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Attn:  Dr.  Leslie  G.  McCracken,  Jr. 
Code  3933A 

Office  of  Naval  Research 
Department  of  the  Navy 
Attn:  Geophysics  Branch,  Code  Ijl6 
Washington  25,  D.  C. 
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Defence  Research  Member 
Canadian  Joint  Staff 
2001  Connecticut  Street 
Washington,  D.  C. 
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Attn:     Dr.  Mervin  J.    Kelly 


Engineering  Library 
University  of  California 
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